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Abstract: The topological interaction arising in interlinked polymeric rings such 
as DNA catenanes is considered. More specifically, the free energy for a pair of 
linked random walk rings is derived where the distance R between two segments 
each of which is part of a different ring is kept constant. The topology conservation 
is imposed by the Gauss invariant. A previous approach (M.Otto, T.A. Vilgis, 
Phys. Rev. Lett. 80, 881 (1998)) to the problem is refined in several ways. It is 
confirmed, that asymptotically, i.e. for large R 3> Rg where Rq is average size 
of single random walk ring, the effective topological interaction (free energy) scales 
(xR A . 
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1 Introduction 



Despite their long history, topologically constrained system such as knots and links 
become an increasingly popular subject of research in statistical physics [U 121 El 
mElEj- Chemically synthesized links (called catenanes) have been considered very 
early pH|H]. In nature, multiply linked rings made of DNA occur in bacteria such as 
Escherichia coli during the replication process as intermediate products 9J. Specific 
topological states (links) can be formed artificially by turning on and off enzymes, 
so-called topoisomerases, which cut and glue the polymeric rings ;9] (see also [HI]). 
These methods have been used to study experimentally the conformational statis- 
tics of specific DNA catenanes [TT], whose elasticity could be studied as well using 
single-molecule techniques [12J. From the viewpoint of polymer physics topological 
interactions due to fixed entanglements remain a challenging problem for systems 
such as rubber networks [T31 114j. 

A fundamental problem in theoretical approaches to linked polymer rings systems 
is the specification of topological constraints. The simplest approach introduced by 
Edwards pQ is a two rings invariant, known as the Gauss integral. It is one of the 
simplest invariants involving a double line integral of a function based on respective 
polymer segment positions of each ring. Thus, it allows to couple explicitly polymer 
conformation and topology conservation. Although being a topological invariant, i.e. 
being independent of the specific polymer conformation as long as the topological 
state (here the mutual linking number) is conserved, it has the disadvantage of 
not being one-to-one. The so-called Whitehead link (being entangled due to self- 
intersections one single ring with itself, see figure GJ and the trivial link (consisting 
of two non-entangled rings) have both the same linking number zero (see jT^j for 
lucid discussion of this point). This result has roused criticism against the use of 
the Gauss invariant in the context of entangled polymer rings, see e.g. |16j . 

The present work responds to these objections in the following way. First, we 
restrict our discussion to random walk rings which are particularly simple to deal 
with within an analytical approach. In this Whitehead link can be disen- 

tangled (justifying a linking number zero), because only the self-entanglement of 
one rings with itself (see figure keeps the link from falling apart. The second 
justification for using the Gauss integral is the fact that the Gauss integral appears 
as the first relevant two-rings invariant appearing within a topological perturbation 
series constructed from averages (vacuum-to-vacuum expectation values, vev's) of 



Figure 1: a) The trivial link: 2 non-entangled rings with linking number 0. b) Two 
entangled rings with linking number —1. c) Whitehead link: two entangled rings 
with linking number 0. The oriented surface bounded by the "thick" ring is pierced 
by the other ring once in a positive sense, once in a negative sense. 
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Wilson loops with respect to the non-abelian Chern-Simons theory |17j . The latter 
has been shown to be equivalent JB] to the HOMFLY knot polynomial which is 
an algebraic invariant usually defined in terms of skein relations, i.e. a recursive 
equation relating different sets of crossings (of the 2D projection of a knot or link) 
giving a recursive algorithm of transforming different knots and links into each other 
|17j . In conclusion, the Gauss integral appears as a quantity to define a minimal 
model for topology conservation which is well defined within in a more complete 
framework just mentioned [TT?] . 

More specifically, the present work revisits the problem of calculating the topo- 
logical interaction for a pair of mutually interlinked polymer rings. This is ac- 
complished by introducing the constraint that the distance between two segments 
each belonging to a different ring be equal to R [201 • Previous work J2U| on the 
problem contained several approximations: the distance constraint was not strictly 
implemented within conformational averages but the distance vector was somehow 
extracted by a pre-averaging approximation. This approximation can be given up, 
as the present work will show. Moreover the discrete nature of the linking num- 
ber is treated within a systematic approach due to Iwata and Kimura [2H] (it has 
recently received some attention in Perhaps not too surprisingly, the anhar- 

monic attraction oc R 4 found in [20] still holds for large distance R/Rg ^ 1 where 
Rg ~ vNP is radius of gyration of a single random walk ring. 

The outline of the paper is as follows. In section 2, we discuss topology con- 
servation using the Gauss invariant together with a distance constraint and define 
conformational averages. In section 3, the (effective) topological interaction between 
chain segments belonging to different rings is calculated: in subsection 3.1, we review 
the topological moment expansion due to Iwata and Kimura [23], in 3.2 we derive 
the topological interaction within the approximation of using the second topological 
moment conditional on a distance constraint M 2 (R), before finally giving an explicit 
expression for the latter in subsection 3.3. Finally we give results in subsection 3.4 
for the topological interaction and discuss its dependence on segment length N . In 
section 4, a conclusion and a brief outlook is given. 

2 Topology conservation with the Gauss invariant 
and a distance constraint 

We consider a very simple system, i.e. a pair of flexible rings which are mutually 
interlinked as in figure [2] and which are subject to an additional distance constraint. 
The latter controls the distance between two segments each of which is part of a 
different ring and is introduced to monitor the effective topological interaction which 
arises due to the topological constraint. 

Based on justifications mentiond in the introduction, the topological constraint 
imposed on the pair of rings - the fact that their are wound around each other as in 
figure [2]- is expressed in terms of the Gauss invariant 



The vector functions r a = r a (s) for a = 1, 2 denote the position in three-dimensional 



$(Ci,C 2 ) =$ 




dr 1 A dr 2 ■ 




(1) 
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Figure 2: Two entangled rings with two segments located on different rings kept at 
a distance R. 



space of segments along the polymer contour, parametrized in terms of s, for each 
ring C a . The Gauss invariant is always an integer for a given topological state re- 
gardless of the deformations imposed on the polymer conformation, i.e. the segment 
position r a (s), which leave the topological state invariant (siehe z.B.JTTj). The in- 
teger which is denoted n here is usually called the linking number and counts how 
often one chain is wound around the other before being closed to a ring (apart from 
a sign which depends on the respective orientation of each ring). 

The problem that we propose to study is the partition sum of a given topological 
state, i.e. a linking number n, for a pair of rings subject to the distance constraint 
as shown in figure El 

Z(n;R) = (5(n,$ 12 )) R (2) 

Inside the average to be defined below a Kronecker delta appears which is equal to 
1, if the arguments are equal, and else. The average is taken w.r.t. to random 
walk ring conformations 

(...) R = jvr\s) Jvr 2 (s)5(r 1 (0)-r 1 (N))5(r 2 (0)-r 2 (N)) 

x S(r\0) -r 2 (0)-R)...e^ 1+ ^ 2 -^- (3) 

Zl2, R 

where 



2P J \ ds f 

and Zi2,n is equal to the numerator on the r.h.s. of Eq.© without the dots. Let 
us note here that actual averages (carried out in appendix A) involve a discretized 
version of the continuous formalism given above. 

The problem Eq.(J2J) has been studied before by Iwata and Kimura [22] and more 
recently by Otto and Vilgis [20] • hi |2"5] . the integrand of the Gauss invariant was 
however approximated by its behavior close to [i -1 — r 2 | ~ 0, an approximation 
which is not necessary as will be shown below. In [20], the Gauss invariant was 
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treated in full generality, but a full implementation of the distance constraint was 
circumvented by a pre-averaging procedure involving the tangent vector densities. 
This approximation can be avoided as will be shown below. The previously obtained 
results jlOj are in fact validated by the subsequent analysis. 

The choice of the random walk ensemble for polymer conformations is motivated 
by the search for an analytical solution. Inclusion of the excluded volume interac- 
tions appears possible in principle (in the context of 0(N) field theories [21]) but 
seems to render the problem of finding an expression for Z(n; R) rather hopeless. 
A derivation of topological interactions from the Gauss invariant treating the latter 
in its most possible generality has not been done up to now even for the simple case 
of random walk conformations (except for a first approach in |20j). 



3 The effective topological interaction 
3.1 Topological moment expansion 

We next proceed to calculate Z(n;R). We rewrite the Kronecker delta inside the 
brackets in Eq.(J2J) by introducing a topological charge, 

Z(n; R) = f p-e^ie***)* = f %^Z{ 9] R) (5) 

J-7T 271_ J-7T ^ 

As already noted in [23], the function Z(g; R) can be expanded in terms of topolog- 
ical moments of order k, ((§n) k )n,- Assuming that Z(n; R) is an even function of 
n, only the even moments k = 2p for p an integer do contribute. This yields 



% r ) = E|^m r ) ( 6 ) 



where 

M 2p (R) = (($ 12 ) 2p )r (7) 

Next we assume that the 2p-correlation can be decomposed into products of p factors 
involving the correlation (($i 2 ) 2 )r so that 



((^ 2 ) 2p )r-^t((^ 2 ) 2 )r (8) 

Z(s;R) = exp(-^M 2 (R)) (9) 

We are thus left with determining the second topological moment given the distance 
constraint R, M 2 (R). Before giving the explicit functional dependence on R , 
we first derive the topological interaction for a given linking number and distance 
constraint in order to show how M 2 (R) appears in the interaction. 



Then one obtains 
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3.2 Topological interaction from the second topological mo- 
ment 



In order to establish the effective topological interactions between segments sepa- 
rated by a distance R — |R| on different rings which are mutually entangled, we 
return to the evaluation of the constrained partition sum Z(n; R) introduced above. 
The second topological moment M2(R) appears - in the approximation Eq.fjHJ) given 
above - inside the conjugate partition sum Z(g;H). It is appears natural now to 
directly perform the integral in Eq.fjSJ). However this leads to negative values for 
Z{n\ R) for even values of n, while odd values of n give positive values, a problem al- 
ready discussed by Iwata and Kimuara[23]. It is traced back there to an unbalanced 
distribution of errors involved in the approximation of breaking up higher moments 
of <&\ 2 as done in Eq.(jHJ). Positive errors are accumulated for odd n and negative 
errors for even n. [2E] 

We therefore proceed like in |25j . using the "continuous Fourier transformation 
(CFT) method", whose basic steps are briefly recalled here. While Z(n; R) is defined 
for discrete n, a continuous function ((t; R) can be obtained by defining 



C(*;R)= -n)Z(n;R) 

n 

Its Fourier transform reads as 

POO 

C{u; R) = / dte' lu \(t- R) 



(10) 



11' 



On the level of Fourier transforms, R) and Z(g; R) are identified for u = g 
which yields 



C(u; R) = Z(u; R) = exp ( — — M 2 (R) 



u 



(12) 



The inverse Fourier transform yields 

*°° du 



C(t; R) 



2tt 



e lu X(u; R) 



2ttM 2 (R 



exp 



t 2 



2M 2 (R) 



(13) 



Now, the constrained partition sum of discrete linking numbers n ^ is obtained 
by integrating 



Z(n; R) 



dt((t; R) 



(14) 



In the present case, one finds 
Z(n;K) = I 



erf 



V^M^R) 



— erf 



n 



V / 2M 2 (R) 



(15) 



where erf(x) is the error function. As noted in the CFT method leads to the 
following problem. Even though £(t, R) is assumed to be a multipeak function in 
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t, the final result just obtained is a smooth function in t. In other words the fine 
structure in t is lost, which is related to the approximations imposed on M 2p (R) for 
large p ^> 1 j2SJ- In the following we will not be interested in the exact quantitative 
dependence of Z(n; R) on n, but rather on the asymptotic behavior for large R = 
|R|. Moreover, the errors contained in the CFT method are evenly distributed 
among different n. 

The effective topological interaction can be derived from Z(n; R) as follows 



(3F(n; R) = -lnZ(n; R) 



(16) 



where (3 = is the inverse thermal energy. If M 2 (R) = M 2 (R) is small and/or 

|n|±- 

n ^> 1 which happens for large R — |R| (see below), such that x± = t _ _ , 2 _ , ^> 1 



then the following expansion will be used 



pF(n; R) = /3F(n; R) 



\n 



2M 2 (R) 



+ In 



n 



^2nM 2 {R) 



0(l/x. 



(17) 



We see that the behavior of the topological interaction (3F(n; R) w.r.t. R, in par- 
ticular for R ^> lyN is dominated by 1/M 2 (R). 



3.3 The second topological moment 

We now proceed to calculate M 2 (R). Using the Fourier integral representation of 
the Gauss invariant 

$12 = ~i I f da [ d5'r;( S )^( S ')e iq - (rl W- r2 ^e^ A % (18) 

VqJo JO Q 

where summation over repeated indices is implied and J . . . — J d 3 q/(27r) 3 . . ., the 
second topological moment reads as 

M 2 (R) = - J jj\ 4 a=1 (£ ds^ e^ x e aTp ^ 

x (r;(si)^( S 3)^(s 2 )^(s4)e iq - (rl(si) - r2(s3)) e ik -( rl ^)- r2 ( S4 ») (19) 

The term in brackets inside the integrals can be rewritten as 

= (ri( Sl )^( S3 )r:( S 2)^( S4 )e^( rl ( Sl ^ r2 ^»e jk -( rl (-)- r2 (-))) R 

= e i ^ +k > R iV(q,k,s 1 , S2 ) J P^(-q,-k,s 3 ,%) (20) 

The factorization on the r.h.s. is due to the factorization of conformational averages 
with respect to polymer 1 and 2 in Eq.(j2J). Details are given in appendix A. We find 
for F Ato .(q,k,si,s 2 ) that 

^(■Sig CT + s 2 k a ) - k a 



F pa (q, k, si,s 2 ) 




7 



X N 
x exp 



i 2 r i 

Sl (q + k) 2 + (s 2 - sjk 2 - — (siq + s 2 k) 2 



6A 

+ [si^s 2 ,q^k] (21) 

The last line on the r.h.s. indicates the operations to be done in order to obtain 
a fully symmetrized expression (see appendix A). Let us note first that the term 

following the factor (j^j on the r.h.s. does not contribute to the second moment, 

because expressions like . . . q^qxe^x ■ ■ ■ vanish. In order to further simplify the 
subsequent analysis, we neglect terms involving 1/N as we assume iV 3> 1 like in 
j2H|- Then F^ a (q, k, s±, s 2 ) reads as 



F^q, k, si,s 2 ) 



— ) $na$(si - s 2 ) exp <{ - g-«i(q + k) 



The second topological moment may then be computed as 



Mo(R) = -2 




/(q + k,iV) : 



q Jk 



q 2 k 2 



i(q+k)-R 



where 



N 



dse-ir su2 



/(u,JV) = 

Transforming k — > u = q + k, one obtains 
M 2 (R) = 



1 — e s 




q- (u-qje iuR / -in u ^ 2 

1 — e e 



u r|u- q 



2 M 4 



f22l 



(23) 



(24) 



(25) 



Introducing u' = \Jl 2 N/ (6)u and X = R/ \Jl 2 N/ (6) and dropping primes, the 
double integral is slightly simplified: 

,iu-X 



M 2 (R) 




q - qj e* 



12 M 4 



1 -e" u 



(26) 



' u r|u-q| 

The integration which involves some approximations affecting the short distance 
behavior (to be discussed below) and which is given in appendix B, yields 



M 2 {R) 



1 



A^- - 7T l F 1 



13 X 2 
2' 2'T" 



-X 2 /4 



(27) 



where = X\Jl 2 N/Q and i_Fi (a; (3; x) is the confluent hypergeometric function. 
The expression involves a cutoff A ~ y/N (see appendix A) whose relevance will be 
discussed below. It is interesting to note what happens when M 2 (R) is integrated 
over all distances between the segments. Then the unconditional second topological 
moment is obtained as 



Mo 



V~ l \ d RM 2 (R) 



V- 1 v^iv76 3 4vr r dXX 2 ^- f A— - it x F x ( --; -; — 
V ' J 8vr 3 V 2 V 2 2 4 

V~H 3 N^ 2 ( Cl N l l 2 - c 2 ) 



-X 2 /4 



(28) 



where in the last line A ~ yiV has been used and where c\ and c 2 are numerical 
factors. Note that the factor 1/V on the r.h.s. of the previous equation is understood 
in the limit V —>■ 00 and N 00 such that p = N/V is a density. This situation 
is given if 2 rings fill out a macrocopic volume. Then the unconditional second 
topological moment has the scaling 

M 2 ~ pi 3 ( Cl N - c 2 X 1/2 ) (29) 

a result, which up to numerical factors C\ and c 2 has been found for the topological 
moment of randomly entangled polymer rings in a dense system where the excluded 
volume effect is ignored |2*2*1 12*1] . Therefore, the cutoff-dependent first term on the 
r.h.s. of Eq. (J27j) appears essential in order to establish consistency with the well 
known result for M 2 for random walk rings to scale to dominant order like oc N. 
We will see below that it affects the effective segment-segment interaction only for 
small up to intermediate deformations of the rings. 

As shown in the previous section, the behavior of 1/M 2 (-R) is particularly im- 
portant to obtain the large R behavior of the topological interaction. Noting that 

^(4I4)=^( 2 '?-t) < 30 » 

and using for large X ^> 1 



one obtains 



4tt\ 1 



7r M^-Aexp(-X 2 /4) + - I 



M 2 (R) V 2 X*J 

2vr 2 X 4 (32) 



for large X = R/^NP/6. It gives the following effective topological interaction 

/3F(n;X)^7r 2 (|n|-0 2 X 4 (33) 

Concerning the X 4 dependence, this result agrees perfectly with the previously 
derived asymptotics for the topological interaction of two concatenated chains based 
on a pre-averaging approximation 



3.4 Results 

Given the expression for the conditional topological moment Eq. (j27j) . we consider 
now the effective topological interaction in more detail. For large X, we use the first 
term in the expansion Eq. (jl7j) . For small X where the Eq. ()17|) is bound to fail, the 
exact expression Eq. (fT5j) for the partition function (based on the second topological 
moment conditional on R or X) is used. In Fig. 01 curves for small X (obtained 
numerically from Eq. (|15jl ) for the topological interaction (3F for linking number 
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n = 1 for various chain lengths N are shown, together with the corresponding 
curves for all X calculated from the first term in Eq.(|17|h whith both sets of curves 
involving Eq. (j27|) for the 2nd topological moment M 2 (R). 




,,c — ■- — — — n 

0.01 0.1 1 10 100 1000 

X 



Figure 3: The topological interactions for linking number n = 1 based on, one, 
the first term in Eq.(|17|) (thick lines), and two, on the full expression Eq.(|15jl for 
the partition function (thin lines), both together with Eq. (j27|) . The thin dotted line 
corresponds to the asymptotic expression Eq.(j3HJ). The inset shows in a double- linear 
plot how the curves based on the full expression Eq. (JTKJ) approach the approximate 
ones (Eq.(|nj)). 

The topological interactions computed from the full expression Eq. (|T5|) smoothly 
approach the approximate ones (valid for large X) based on Eq.(|17|). Whereas for 
large X the topological interaction is independent of the chain length TV approach- 
ing the asymptotic form Eq. (j3l?]) . such a dependence appears for intermediate and 
small X. To be precise, the interaction becomes weaker for larger N for a given in- 
termediate X, just below X = 10. This effect might be due to the discrete polymer 
model used here. Let us consider two systems of two concatenated rings which are 
similar in their conformation but which have different chain length Ni and N 2 with 
N 2 > N\. Then the system with chain length N 2 suffers a smaller entropy loss due 
to the topological constraint. Thus the topological interaction is reduced. 

However, for very small X the topological interaction does not show a monoton- 
uous dependence on N but passes through a minimum as N is increased from 10 3 to 
10 9 . Moreover, for very large N a minimum of the topological interaction as function 
of X starts to develop. A strong repulsion of segments belonging to different rings 
at small X found in [20] is not confirmed by the present calculation. On the other 
hand, the present work depends on approximations (see appendix B) which possi- 
bly influence the short distance behavior of the topological interaction. Therefore, 
no conclusive statement is claimed for the topological interactions derived from the 
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Gauss invariant on small scales X. A more complete analysis of the short- distance 
behavior is necessary and the subject of future work. 

For extreme elongations X ^> 1, one might expect for the topological interaction 
a crossover from X 4 to X 2 as claimed in [21], based on the argument that when 
two random walk rings linked with linking number 1 are fully stretched, they form 
an effective random walk chain. Apparently, this limit cannot be reached by the 
present approach. The argument just given seems to assume, however, some local 
gluing and some partial inextensibility of the chains which then leads to a mutual 
attachment of the rings due to the topological constraint. Each ring conformation is 
however treated here w.r.t. to the ensemble of random walks. Therefore the effective 
conformation of a single random walk chain is not reached. 

4 Discussion 

In the present work, the topological interaction between segments belonging to dif- 
ferent polymer rings which are mutually interlinked is considered. Two segments are 
singled out and are kept at a distance R, while the conformations of each ring is av- 
eraged over w.r.t. a random walk ensemble. Supplementing a previous calculation 
20J, the present analysis introduces several refinements. First, the distance con- 
straint is implemented from the very beginning. The discrete nature of the linking 
number is treated systematically following 25\. 

The topological interaction behaves asymptotically, i.e. for X ^> 1 where 
X = R/y/NP/6, as ~ X 4 as found previously [2D]- For smaller X < 10, the 
topological interaction depends on the chain length N. For intermediate distances 
X, an argument related to the reduced entropy loss due to the topology conserva- 
tion for increasing N can be given. For small X, no final conclusive statement on 
the topological interaction can be given, as the approximations used for the confor- 
mational averages in order to obtain analytical results possibly affect small length 
scales. We remark however that the present result is consistent with calculations of 
the unconditional second topological moment [221 121] • 

The present work may be considered as a further step forward to understand 
the elasticity of DNA catenanes (interlinked rings) which can be studied in principle 
using single-molecule techniques [12]. A first step was the previous study [20] in 
order to understand the conformational statistics of DNA catenanes ^T] . The simple 
random walk model used here is of course not realistic in this context. Apart from 
the excluded volume problem, DNA molecules are semiflexible chains [2Z] which 
requires the Kratky-Porod model [2H] which however leads to analytical problems. 
Nonetheless, for long chains, self-avoiding walk behavior is recovered. 

Possible future directions of the problem presented here not only concern to use 
more realistic conformational models relevant to the study of biomolecules. Future 
work should aim at a more accurate treatment of conformational averages, in par- 
ticular the inclusion of terms which were omitted here for large N (which amount to 
neglect the closure of chains). Another important problem is the issue of localiza- 
tion: i.e. the hypothesis, that a number of segments are localized close to a region 
containing the link, motivating concepts such as slip-links (see e.g.|29j) which have 
found renewed interest [3U] . 



11 



Acknowledgements: Financial support by the DFG under grant Zi209/6-l is 
gratefully acknowledged. 



Appendix A: The correlation function 7^ 



Following [2H], we calculate the following correlation function of Eq . (|20p 



1 fMUO 



r 1 l' e „V 2 l'Q.^ i( l-( rl (si)-r 2 (s3))^k-(r 1 (s2)-r 2 ( S 4)) 



R 



(34) 



To do so, we introduce a discretized random walk r 1 (s) = Rj + J2i=i D ; sucn that 
the measure of the path integral in Eq.fjSJ) becomes 



Vr\s)^ [d*Rlf[[d 3 b} 

J i=i 



(35) 



Then the correlation function given above is derived from the following generating 
functional 



G ({h]}, {hf }; q, k; a u s 2 , s 3 , s 4 ) = ( e ^^^e A <^-^ 



{R},{hJ},{h?} 



(36) 



The average is understood as 

G ({h}},{h^} ] c l XsuS2,s 3} s 4 ) = ^- J d 3 Rll[ J d 3 b}J d 3 RlH J d% 2 



X 



^(E b ') 5 (E b ^( R o- R o- R ) ex p 

i i 

x exp [-i(q + k) -R* +i(q + k) -R^] 



a=l,2 i 



x exp 



*2 



«4 



i=l 



1=1 



1=1 



8=1 



(37) 



Products and sums in the equation above marked by an index i only are supposed 
to cover the range i = 1..N. The integrations w.r.t. to Rq for a = 1,2 in the 
numerator and denominator on the r.h.s. of the equation above are easily carried 
out to give 

G({hi},{h.2};q,k; Sl , S2 , S3 ,S4) = e l ^ R F ({h]},q,k, Sl , s 2 ) 

x F({h 2 },-q,-k,s 3 ,s 4 ) (38) 



where 
F({hi},q,k,si,s 2 ) 



Z 

x exp 



2l< 



$>,) 2 + ]>> • b, 



81 



S-2 



iq ■ ^ bi - ?k • ^ hi 



i=l 



1=1 



(39) 
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The partition sum Z\ is the single ring partition function devided by a volume factor 
(originating from the integrations w.r.t. to Rq for a = 1, 2 which have already been 
carried out). Now one has to choose an order for si, s 2 which are now considered as 
being discrete out of 1..N, say si < s 2 . Then the integrations are easily carried out 
to give 



F({h;},q,k, s 1 ,s 2 ) 



exp 



— ( £> + Sl q+a 2 k 



*2 



2V 



- nr (h l + q +k) 2 + £ (h ? , + k) 2 + £ h 2 



i=si+l 



i=s 2 +l 



(40) 



The case s\ > s 2 is obtained by interchanging s\ <-> s 2 and q <-» k. Now, the 
correlation function 7 AU/OT is obtained as follows 



1 



d d d d 



-i) 4 dhj^ dhl >S3 dhl >S2 dhl Si 



G ({ h !}> { h i }; q> k ; s i, s 2, s 3 , s 4 ) | h i=o,h?=o 

(41) 

Using Eq.(|55jl together with Eq. (}4T)]) one obtains Eq.(|2Djl with Eq. (j2*TJl for j^ vaT given 
in the main text. 



Appendix B: Approximate evaluation of integrals 

Defining the integral w.r.t. q in Eq.([26|) as /(u) one obtains 
where 

W = /(U) = / ! 2 |1 U " ' ft = WY2 B ^l^ 3/2) [ A dq ff-f (43) 
y q r|u-qr (2tt) 2 J (g 2 + u 2 ) 2 

Note that the integral in the last eq. needs to be cut off by a parameter A which is 
related to the cutoff in momentum space A' via A = A'^PN/G and thus proportional 
to VN if A' ~ Evaluating the g-integral yields for l(u) 

1 71 ( 71 Au 2 \ 

l{u) = j^y 2 - i^-u -2A + 3narctan(A/ M ) - -^-^ j (44) 

/ 2\ 2 2 

In order to make analytical progress, let us approximate ( 1 — e~ u J /u 3 ~ we~" in 

Eq . (J26j) . The error involved in this approximation concerns large u, and the short 
distance behavior of M 2 , which is discussed in the main text. Then the second 
topological moment is given by the expression 

16 1 z* 00 

M 2 (R) = -j— — / dul(u)usm(uX)e- u2 (45) 
(2vr) 2 X J 
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Inside the integral l(u) is approximated, for A ^> 1, as 

'W= ! ^pK" aA+a ™"S (46) 

where arctan(A/u) ~ tt/2 for A ^> 1. Using the approximation above, evaluation of 
the integral w.r.t. u involves the following integrals. The term of order in u gives 
the integral 

h = T duu sin{uX)e- u2 = ^ e -* 2 /* (47) 
Jo 4 

The term linear in u in l{u) yields the integral 

h = jf ' duu"sm(uX)e- u2 = fe^^ (-1; |; (48) 

while the term of order 2 in u involves evaluating 

/°° X ( S A^ 2 \ 

^ M 3 sin( M X)e- 2 = -e- x2 / 4 r(5/2) ^ ( -1; ~; —j (49) 

The function iF\ f — 1; |i ^1 explicitly gives 1 — X 2 /6. However, we will neglect 

the u 2 term in l(u) because it is multiplied by a factor 1/A (see Eq.(jlHj)). Collecting 
terms, one obtains for the second moment 

where R = Xy/PN/6 and 1F1 (a; ft; x) is the confluent hypergeometric function (see 
Eq. (|27j) in the main text). 
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